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Abstract 

We consider a lattice discretization of a covariantly gauge-fixed abelian gauge 
theory. The gauge fixing is part of the action defining the theory, and we study the 
phase diagram in detail As there is no BRST symmetry on the lattice, counterterms 
are needed, and we construct those explicitly. We show that the proper adjustment 
of these counterterms drives the theory to a new type of phase transition, at which 
we recover a continuum theory of (free) photons. We present both numerical and 
(one-loop) perturbative results, and show that they are in good agreement near this 
phase transition. Since perturbation theory plays an important role, it is important 
to choose a discretization of the gauge-fixing action such that lattice perturbation 
theory is valid. Indeed, we find numerical evidence that lattice actions not satisfying 
this requirement do not lead to the desired continuum limit. 

While we do not consider fermions here, we argue that our results, in combination 
with previous work, provide very strong evidence that this new phase transition can 
be used to define abelian lattice chiral gauge theories. 
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1 Introduction 



In this paper we continue our investigation of the gauge-fixing approach to the con- 
struction of lattice chiral gauge theories. In this approach, gauge invariance is broken 
both through the gauge-fixing terms and through the fermions. This requires adding a 
complete set of counterterms to the theory, in addition to the gauge-fixing terms, and 
these counterterms will need to be tuned. Showing that this can be done corresponds 
to demonstrating that the phase diagram contains a continuous phase transition which 
can be employed to construct the desired continuum chiral gauge theory [jl], ||. 

Here, we will restrict ourselves to the abelian case. This avoids many of the subtle 
questions concerning Gribov copies which arise in the non-abelian case. In particular, 
it makes it possible to drop the ghost sector from consideration [^, while still testing 
many of the key elements in this approach to lattice chiral gauge theories. 

We will employ (a generalization of) the lattice gauge-fixing action proposed in 
Ref. 0]. Since we wish to maintain close contact with standard weak-coupling per- 
turbation theory, we consider a lattice version of the Lorentz gauge. (Other gauges 
may work as well, but we believe that it is important to restrict oneself to a renor- 
malizable gauge.) On the lattice, i.e. in the regulated theory, there is no gauge or 
BRST symmetry, and both the transverse gauge fields and the fermions will couple 
to the gauge degrees of freedom, represented by the longitudinal part of the gauge 
field. 

This leads to the following simple questions, both of which can be addressed 
without simulating the full theory including all its dynamical degrees of freedom: 1) 
when we turn off the transverse part of the gauge field, do we obtain a theory of free 
(chiral) fermions in the correct representation of the gauge group (in the abelian case, 
with the correct charges), decoupled from the longitudinal modes? and 2) when we 
turn off the fermions, do we obtain a theory of free photons, again decoupled from 
the longitudinal degrees of freedom? 

It is well known that (most) small perturbations of the gauge-invariant compact 
lattice formulation of a U(l) gauge theory do not change the nature of its (weak- 
coupling) continuum limit (they correspond to irrelevant directions) ||^. However, 
previous work has shown that it is generically not possible to construct chiral gauge 
couphngs near such a continuum limit (for reviews, see Refs. f^, 0])- In our approach, 
gauge fixing plays an essential role in the construction of the theory. This means that 
the coupling in front of the gauge-fixing action has to be large enough to bring us 
to a new type of continuous phase transition in the phase diagram, at which both 
questions above can be answered affirmatively. 



We have addressed question 1) in previous work |I0|, |TT], [T^. We showed 

that, in the "reduced" model, in which only the fermions and the longitudinal gauge 
fields are kept, indeed a new type of phase transition occurs. At this phase transition, 
a continuum limit can be defined in which free chiral fermions with the correct charges 
emerge, decoupled from these longitudinal degrees of freedom. Here, we address the 
second question. We turn off the fermions, and we ask whether this new phase tran- 
sition survives when the transverse degrees of freedom are present, and in particular, 
whether it allows us to construct a theory of free photons at the transition. We would 
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like to emphasize that, even though, without fermions, this is not a new continuum 
theory, our critical point corresponds to a new type of universality class. This is the 
key element that allows us to couple fermions chirally to the gauge fields, without 
running into the problems which made many previous attempts unsuccessful. 

In addition to this fundamental question, we also address a more technical issue. 
It was argued in Ref. that one has to be very careful with the precise definition 
of the lattice version of the gauge-fixing action. A naive "standard" discretization of 
the Lorentz gauge- fixing action (1/2^) / d^x {d ■ A)"^ will lead to the occurrence of a 
dense set of lattice Gribov copies (with no continuum counterpart). This corresponds 
to a large class of uncontrolled, rough fluctuations in the lattice theory, and may well 
spoil the existence of the critical point we are after. The lattice gauge-fixing action 
proposed in Ref. [§] does not suffer from this problem, and it is this action that we 
have used in our previous work. Here, we introduce a one-parameter class of lattice 
gauge-fixing actions, which interpolates between the naive discretization and the one 
of Ref. 0]. This corresponds to adding a direction to the phase diagram, and we 
explore the dependence of the phase structure on this new direction. 

The organization of this paper is as follows. In Sect. 2 we give the full action for 
a gauge-fixed U(l) gauge theory, including a complete set of counterterms. In partic- 
ular, we introduce the parameter f, which interpolates between the naive gauge-fxing 
action and that of Ref. 0, and we discuss the above mentioned lattice Gribov copies. 
We argue that standard lattice perturbation theory should be valid as long as r is 
large enough. In Sect. 3 we present analytic results in preparation of a high-statistics 
numerical study of this model. We first explain the nature of the new phase transition 
(which we will denote as the "FM-FMD transition") from the classical potential, and 
then provide a simple-minded mean-field analysis of the model. Since there is good 
qualitative agreement between mean field and our numerical results, we also give an 
overview of the structure of the phase diagram at this stage. We end this section with 
a calculation of the one-loop lattice photon propagator, and use it to determine some 
of the counterterms at one loop. We also consider a (composite) scalar two-point 
function. Then, in Sect. 4, we present our numerical results, which constitute the 
main part. First we discuss in detail how we determined the phase diagram. After 
that, we compute vector and scalar two-point functions numerically, and compare 
them with perturbation theory in order to determine whether we do indeed obtain a 
theory of free photons at the FM-FMD transition. Finally, we summarize and dis- 
cuss our findings in the last section. There are two appendices, containing various 
technical details. 

We would like to end this section with mentioning that, recently, a gauge-invariant 
construction of (anomaly-free) abelian chiral gauge theories on the lattice has been 
proposed, based on a Dirac operator satisfying the Ginsparg- Wilson relation |13[ (for 



the non-abelian case, see Ref. The problems with the violation of gauge invari- 



ance discussed above do not apply in this case (they might if a gauge non-invariant 
approximation of the Dirac operator is used, however). An essential ingredient is 
that the fermion measure includes a gauge-field dependent phase factor which is de- 
termined by requiring the theory to be gauge invariant and local. Sofar, however, no 
explicit expression of the fermion measure was given, making this approach as yet 
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not suitable for a numerical investigation. 



2 The Model 

The central idea of the gauge-fixing approach is to control the gauge degrees of free- 
dom by a gauge-fixing procedure. The starting point is the gauge-fixed action in the 
continuum, the "target theory." Correlation functions of the target theory in the con- 
tinuum satisfy Slavnov-Taylor identities, as a consequence of the gauge symmetry. 
For an abelian gauge group, which will be the subject of this paper, the target action 
in the continuum is of the form, 

= S,^g{A^) + S,,f{A^; ^l, i^R) + ^c,g.f.(^M) . (2-1) 

where Sc,q{A^) = \ j d^x F"^^ designates the gauge action, Sc^F^A^^ipL^ipR) is the 
chiral fermion action, and iS^gi . (^m) gauge- fixing action. Here, we will consider 
the Lorentz gauge, which is renormalizable, and therefore allows us to study the 
(relevant part of the) phase diagram in perturbation theory. No ghosts are needed 
in the abelian case, and henceforth we will not introduce any ghosts on the lattice 
either ||^. We then have 

Sc,,.i{A,) = ^ld'x(Y.d,A,) ' (2-2) 

where ^ is the gauge-fixing parameter. The goal is now to transcribe the target theory 
defined by the action (|2.1|) to the lattice using compact lattice link variables and the 
Haar measure as integration measure in the path integral. 

The fermion action leads to violations of the Slavnov-Taylor identities if we use 
any one of the standard lattice fermion formulations, such as Wilson, staggered or 
domain-wall fermions, which are in conflict with chiral gauge invariance. Moreover, 
we will consider a class of gauge-fixing terms, dependent on a continuous parame- 
ter r. They will lead to additional violations of the Slavnov-Taylor identities. The 
Slavnov-Taylor identities are restored in the continuum limit by adding a finite num- 
ber of counterterms to the action, and appropriately tuning their coefficients in 
the continuum limit. 

In this paper, we will drop the fermion sector and, as mentioned in Sect. |I|, fo- 
cus on the question whether the lattice discretization of the target action with U(l) 
symmetry provides a valid formulation of a theory of free photons on the lattice. The 
lattice action is formulated in terms of the compact link variables Ufj_x = exp{iagAf^j:), 
with g the gauge coupling and a the lattice spacing (which we will set to one through- 
out this paper). The lattice action is then given by the expression 

S = Sg{U) + S,,aU) + S,,,XU) . (2.3) 

where Sg{U) is the gauge action, S'g.f.(t/) the gauge fixing, and SctX^) the countert- 
erm action on the lattice. For the lattice transcription of the gauge action Sc^qIA^j^) 
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Figure 1: Example of a lattice Grihov copy. Shown is only a two-dimensional projec- 
tion of the four- dimensional lattice. The U^x fields on all the links attached to x = xq 
(fat lines) are equal to —1 and all other U^x fields are equal to one. 

we employed the standard plaquette action, 

Sg{U) = E E {1 - Re U,^x} , (2.4) 

where Ufj,i,x = Ufj^xU^x+jiUj^x+pUlx the usual lattice plaquette variable. The lattice 
transcription of the gauge-fixing action is more subtle. A naive lattice discretization 
of Eq. (O) leads to 



X \ II ) 

where 

V^x = Im U^x , (2.6) 

^ - 2^ • P-^' 

and A^V^a; = V^x — Vfix-p. (A^ is the backward nearest-neighbor lattice derivative). 
The problem with the naive lattice discretization is that the classical vacuum of the 
action SdU) + S0"%U) is not unique §. It is easy to see that SdU) + S0"%U) 
has absolute minima for a dense set of lattice Gribov copies Unx = Qx ^ ol+fij of the 
classical vacuum Uf^x = 1 for particular sets of Qx- An example of such lattice Gribov 
copies is displayed in Fig. |l|, where we set ^f^, = — 1 at site x = Xq and gx = +1 ai all 
sites X ^ Xq. The Ufj^x fields on the links attached to the site Xq (fat lines) are equal 
to —1 whereas the U^x fields on all other links (thin lines) are equal to +1. These 
lattice Gribov copies are a high momentum lattice artifact with no counterpart in the 
continuum. (They should not be confused with continuum Gribov copies, which are 
a long-distance phenomenon.) It is clear that perturbation theory around only one of 
the absolute minima of Sg{U) + S'^f^'^{U), in particular the classical vacuum U^x = 1, 
may not give a valid description of the theory, and the phase diagram could (and will) 
turn out very different from what one would expect from naive perturbation theory. 

It is however possible to remove the unwanted lattice Gribov copies of the classical 
vacuum (and therefore, by continuity, of field configurations perturbatively close to 
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the vacuum) by adding a higher dimensional operator to the gauge-fixing action in 
Eq. ( |2.5| ). This procedure is similar to Wilson's idea of removing the species doublers 
of the naive lattice fermion action by adding an operator of dimension larger than 
four (the Wilson term). Such higher dimensional operators do not affect the small- 
momentum behavior of the theory, but can be used to change the behavior at large 
momenta. 

The gauge-fixing action we will use in this paper is given by the expression 

S,,XU) = S^J:^^{U)+rK {Cx + Cif-Bl^ , (2.8) 



where 



a = ^n,,(f/), (2.9) 

y 

B. = Y. r"o^^'^ )' > (2.10) 



and 

^{UU = E {U,x 5.,+f.,y + Ul-i. 5x-A,. - 2 Ky) (2.11) 



is the covariant lattice Laplacian. In Eq. ( |2.8| ) we have multiplied the higher di- 
mensional operator by a new parameter r which can be viewed as the analog of the 
Wilson parameter r that multiplies the Wilson term. It can be shown that the action 
Sg{U) +Sg,i.iU) has, for f > 0, a unique absolute minimum at U^x = 1 so that, for 
r > 0, standard perturbation theory in g is valid. The gauge-fixing action provides a 
continuous interpolation between the naive gauge-fixing action ( |2.5| ) (f = 0) and the 
gauge-fixing action at r = 1 which was introduced previously in Ref. @] and was used 
inRefs. [|,|,|10|,|n],|12|, |l|]. 

We will study the r dependence of the phase diagram of the purely bosonic 
theory, and explore the effects of the lattice Gribov copies on the phase structure. 
Obviously, lattice Gribov copies introduce rough gauge degrees of freedom. It is 
therefore conceivable that lattice Gribov copies at small r give rise to a disordered 
or symmetric phase. It has been argued earlier that a chiral gauge theory cannot be 
obtained in such a phase |l^. This would teach us that, in the gauge-fixing approach 



to lattice chiral gauge theories, one needs to choose f of order one. 

Finally, we have to specify the counterterm action. It contains all relevant and 
marginal operators which are allowed by the exact lattice symmetries . In the case 
of the U(l) gauge theory there are six such terms: 



fix X flU 

- ^2 EE Kim u,,)' - A3 E (e u^x) 

X fl X \ fl J 

- A4 E fE(Im - A5 EE (Ini U,xf . (2.12) 

X \ fl / X fl 
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The term proportional to k is a mass counterterm for the gauge field. It is the 
only dimension-two counterterm. All other counterterms are of dimension four. An 
expansion in g shows that the counterterms with coefficients Ai, A2 and A3 are wave- 
function renormalization counterterms. The terms with the coefficients A4 and A5 
are, to leading order in quartic in the gauge potential, and are needed to eliminate 
photon self-interactions. The coefficients of the counterterms have to be tuned such 
that the Slavnov-Taylor identities are satisfied in the continuum limit (c/. Sect. |3.3|) . 
Note that terms such as X)z T.pii.^^^^ U^x) J^uO-^ U^xY do not occur because of the 
symmetry f/^^. —>■ f/^^.. However, such counterterms will have to be considered once 
the theory is coupled to fermions. 
The path integral is defined by 

Z = J DU exp{-S{U)) (2.13) 
= I DUD<P exp{-Si<PlU,,<Px+^)) , (2.14) 



where the action S{U) is given in Eq. ( |2.3|) , and the integration over the link variables 
is performed with the usual Haar measure. Only the gauge action in Eq. (|2.3| ) is 
invariant under the local gauge transformation 

U,. - U^,x = 9.U,xgU^ . (2.15) 

In Eq. ( p.l4|) we made the integration over the longitudinal gauge degrees of 
freedom (p^ explicit. It is easy to see that Eq. ( |2.13| ) follows from Eq. ( p.l4[ ) by 
performing the gauge transformation ( |2.15| ) with g^ = (pl and using the fact that 



/ D(j) =10]. The longitudinal modes can be viewed as group- valued Higgs fields (p^- 
We will therefore denote in the following S{4>lUfj,x4>x+fi) as the action in the "Higgs 
picture" and S{U) as the action in the "vector picture." Both pictures are equivalent 



and every observable in the vector picture has, according to Eq. (|2.15|) , a counterpart 
in the Higgs picture. The path integral in the Higgs picture is invariant under the 
local transformation 

4>x hx(j)x , U^x hxU^xhl^f, . (2.16) 

Finally, we mention here that in d space-time dimensions the path integrals ( p.l3|) 
and ( |2.14 ) are invariant under the additional discrete symmetry (with all Aj = 0) 



Uf,x ^ -Uf,x , K^-K-Sdnr, k^k, r^r, g g , (2.17) 

which implies that we can restrict ourselves to the k > —Adnr region of the {k, k, r) 
phase diagram. In the following we will refer to 

K = ~4dKr (2.18) 

as the symmetry surface, or symmetry line (point) if k or (and) f are kept fixed. 

At R = and Aj = 0, i = 1, . . . , 5, the model reduces to the U(l) gauge-Higgs 
model with group-valued Higgs fields. The resulting {k, g) phase diagram corresponds 
to the K = hyperplane of the four-dimensional (k, g, k, r) phase diagram (for H = 
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the model is independent of f). In the hmit g ^ the phase diagram of the gauge- 
Higgs model reduces to that of the four-dimensional XY model. This phase diagram 
of the XY model contains three phases: a ferromagnetic (FM) or broken phase at 
K > kfm-pm ~ 0.15, a paramagnetic (PM) or symmetric phase at — kfm-pm < «^ < 
/^FM-PM, and an antiferromagnetic (AM) phase at k < — kfm-pm- The k > region 
is mapped by the symmetry (|2.17|) into the k < region. This implies that the FM 
phase is mapped into the AM phase and vice versa. The FM-PM transition, and, 
because of the symmetry (|2.17| ), also the PM-AM phase transition are second-order 
phase transitions. At g > the FM phase turns into a Higgs and the PM phase 
into a Coulomb phase, which we will still denote as FM and PM phases, respectively. 
There is strong evidence that also the FM-PM phase transition at > is of second 



order IT^. The spectrum in the FM phase contains a massive vector boson with 
quantum numbers J^^ = 1 and a massive Higgs boson with quantum numbers 
J^*^ = 0^+. The masses of both particles are expected to scale when k is tuned 
towards the second-order FM-PM phase transition. The spectrum in the Coulomb 
phase contains only a massless photon if one keeps away from the FM-PM phase 
transition [T^. The model in the Coulomb phase provides a valid formulation of free 



photons on the lattice but, as we explained before, because of the strongly fluctuating 
longitudinal gauge modes, it is not possible to formulate a chiral gauge theory with 
a gauge non-invariant fermion action in that phase. The spectrum near the FM- 
PM phase transition contains, apart from the massless photon, also positronium-like 



bound states of scalars |T^. At large g and small n the phase diagram contains a 
confining (CF) phase, which is separated from the PM phase by a phase transition, 
located at (7 ~ 1. 



In Refs. p, 1^, [ly, |TT|, |T^, we have studied the U(l) model with Aj = 0, i = 1, 



in the reduced limit where the gauge coupling is tuned to zero, while keeping k = 
l/{2g'^^) fixed (reduced model). In this limit the transverse gauge degrees of freedom 
are gone and only the dynamics of the longitudinal gauge degrees of freedom remains. 
The path integral of this reduced model is given by Eq. ( |2.14| ) with all link variables 
set equal to one. The phase diagram of the reduced model at r = 1 contains, apart 
from the FM and PM phases, also a new type ferromagnetic directional (FMD) phase 
where the vector field Im0j,0a;+^ condenses and hypercubic rotation invariance on the 
lattice is broken (c/. the f = 1 plane of Fig. |^, which is very similar to the f = 1 
reduced-model phase diagram). 

The phase transition between the FM and FMD phases, which, as we will see, is 
also present in the full U(l) model, plays a crucial role in the gauge-fixing approach. 
While the photon is massive in the FM phase, one obtains a massless photon by tuning 
K, from within the FM phase, towards this FM-FMD phase transition. Tuning all 
other counterterms as well should then lead to a theory of free, relativistic photons. 
We should note here that k is the only relevant counterterm parameter, whereas 
all the Aj are marginal. We therefore may expect that only /t needs to be tuned 
nonperturbatively, while, with a given numerical precision, a one-loop, or even tree- 
level determination of the Aj will suffice. Indeed, the results which we will present in 
the following sections provide evidence of this. 

With chiral fermions, we demonstrated that, in the reduced limit, the fermion 
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spectrum contains only the desired chiral states |T^, |TT], |T2|. 

The phase diagram at small g is expected to look similar to the phase diagram of 
the reduced model, since the transverse components of the gauge fields are still very 
small. We will show in the following sections that even at a relatively large value of 
g {g = 0.6) the phase diagram is qualitatively very similar to the one for the reduced 
model, if f f« 1. 



3 Analytical results 

In this section we present our analytical results. Subsect. PTTI deals with the constant- 



field approximation, which already gives some insight into the phase structure. In 
Subsect. p.2| we will determine the {k, k, f) phase diagram at fixed small g in the mean- 
field approximation. The counterterm coefficients k, Ai, A2 and A3 are calculated to 
one-loop order in perturbation theory in Subsect. p.3| . In Subsect. p.4| we analyze the 
spectrum near the FM-FMD phase transition in perturbation theory. 

3.1 Constant-field approximation 

In the constant-field approximation we set 

U^^ = exp{igA^) , (3.1) 

where is a space-time- independent vector potential. After inserting ( |3.1| ) into 
Eq. ( p.3|) all terms which contain derivatives of the gauge field vanish, and we obtain 
an expression for the classical potential density. Expanding the resulting expression 
in powers of g, we find 



Ki(A^) 



(3.2) 



where ". . ." represents terms which are of higher order in g^. A massless non- 
interacting photon is obtained for k = 0, A4 = and A5 = 0. The relation k = 
fi^FM-FMD = defines a critical surface in the three-dimensional (k, k,, f) phase dia- 
gram where the photon mass vanishes. We will see later that the critical coupling 
/tFM-FMD('5; ^) is shifted away from zero by perturbative corrections. 

The minimization of the classical potential density (|3.2|) for A4 = A5 = shows 

that 

{g A^) = , for K > kfm-fmd , „x 

{g Af,) = ±{\k- /«fm-fmd|/(6 k r)f^'^ , for k < kfm-fmd , 

for /i = 1,...,4 (see Ref. [Q, which also deals with the case A4 7^ 0, A5 7^ 0). 
This implies that k = kfm-fmd = corresponds to a phase transition between the 
FM phase where (Af^) vanishes and the gauge boson has a nonzero mass, and the 
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FMD phase with a nonvanishing vector condensate (A^). The hypercubic rotation 
invariance on the lattice is broken in the FMD phase by the nonvanishing vector 
condensate. 

The constant-field approximation is supposed to provide a satisfactory description 
of the model only at large /5 ~ 1/(7^, where strongly fiuctuating gauge configurations 
are suppressed by a small Boltzmann weight and the smoother gauge configurations 
are well approximated by a constant field. Indeed, as we will see, a perturbative 
expansion starting from the constant-field approximation can be developed by ex- 
panding = ^wi'^gA^x) and setting R = l/{2^g^), with fixed. The situation 
is different at small R where rough longitudinal fields are not sufficiently suppressed. 
This picture is confirmed by our previous investigation of the reduced model at f = 1, 
where we find indeed a FM-FMD phase transition at large k (consistent with the con- 
stant field-approximation), but an FM-PM phase transition at small k. The emer- 
gence of the disordered PM phase at small R is due to the dominance of rough gauge 
field configurations. Equation ( p.2|) shows that the classical potential density depends 
only on the product Rr. This suggests that a small value of f has the same effect as 
a small value of R and that the "small k" -region around f = might indeed be filled 
with a PM phase. This conjecture will be confirmed by our mean-field analysis and 
Monte Carlo simulations. 



3.2 Mean Field Calculation of the Phase diagram 

In the following we set again Aj = 0,i = l,...,5 and determine the phase boundaries 
in the {K,R,r) space. The mean-field analysis for gauge theories on the lattice is 



ambiguous as it is in confiict with local gauge invariance [|T^ . For gauge-Higgs models 
the mean-field calculation in the vector picture leads to a wrong phase structure at 
small values of k, whereas the mean-field approximation in the Higgs picture leads to 
a phase diagram which, at least qualitatively, complies at small k with Monte Carlo 
simulations. Since we are interested in the phase diagram at small k, we performed 
our mean-field analysis in the Higgs picture. Here, we will only describe the results, 
while relegating all the technicalities to Appendix A. 

All the results we are going to present in the following were obtained at g = 0.6, 
scanning the three-dimensional {k, R, r) phase diagram at g = 0.6 by keeping either 
R or r fixed. (We also performed a few scans of the four-dimensional {g,K,R,r) 
parameter space in the g direction at f = 1, k = and several values of R and 
find that the transition to the confining phase occurs always at values of g which 
are larger than 0.6.) The resulting two-dimensional sections through the (k, k, f)- 
phase diagram are displayed in the left columns of Figs. Q-^. The corresponding 
phase diagrams obtained from the Monte Carlo simulations are displayed in the right 
columns. They will be discussed in more detail in Sect. but a first glimpse shows 
that they are, in most cases, similar to the mean-field phase diagrams. 



The dash-dotted lines in Figs, g-]^ mark the symmetry line (|2.18|) ; we only de- 
termined the phase diagram above that line. The error bars in the mean-field phase 
diagrams mark the distance between two successive points in our scans of the phase 
diagram (c/. Appendix A). 
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Figs. I^a and b show the [k, k) phase diagram at f = 1. There is a PM phase 
for small \k\ and small \k\, and for large k there is a phase-transition line separating 
an FM and an FMD phase. The plots in Figs. show that the situation remains 
qualitatively the same when r is lowered. The phase diagram at f = contains at 
small K and k a PM phase, as for f > 0. The mean- field calculation, however, predicts 
that, at f = 0, the region at large k and small \k\ is filled by an FMD phase, which is in 
conflict with our Monte Carlo simulations (plot in the right column), which show that 
the PM phase extends to very large k with no sign of an FMD phase. We therefore 
believe that the FMD phase at large k is an artifact of the mean-field approximation 
(which tends to favor ordered over disordered phases). Figures show three (/t,r) 
phase-diagram plots at fixed k, values. The phase diagrams at k = 0.05 and 0.2 look 
very similar to those at fixed f > 0. For fi;<0<forf<0<Kwe find an FM-AM 
phase transition which coincides with the symmetry line, Eq. ( ^.181 ). Figure ^ shows 
again that the mean-field calculation at /5 = 0.8 does not lead to a PM phase at small 
r. 

We have compiled the Monte Carlo results for the {k, k, r) phase diagram into a 
schematic graph, shown in Fig. |[ The FMD phase at large k and f = is not shown 
in that graph since, as we said, there is no evidence from Monte Carlo simulations 
that the FMD extends down to r = 0. Fig. | shows that the FM-PM and FM-FMD 
phase-transition sheets are separated by a tricritical line where three phases (FM, PM 
and FMD) meet; we will call this the FM-FMD-PM line. Similarly, also the FM-PM 
and FM-AM phase transition sheets are separated by an FM-AM-PM tricritical line 
(not shown in Fig. |^). The projections of these tricritical lines onto a constant-K 
plane are shown in Fig. |^ (mean field) and in Fig. 0b (Monte Carlo). We see that 
the discrepancy between the mean field and Monte Carlo results at f = correlates 
with the fact that the mean-field and Monte Carlo locations of these tricritical lines 
differ slightly. In mean field, the PM phase ends at k ^ 0.3, and the FM-FMD-PM 
and FM-AM-PM lines merge into an FM-FMD-AM tricritical line, whose projection 
approaches the r = axis for k = oo. 

In the mean-field approximation we find that the FM-PM transition is second 
order at /5 = (this is the "standard" gauge-Higgs model). It is still second order at 
K > and r = 1, but changes into a first order phase transition when f is lowered. 
Similarly we find that the FM-FMD transition is of second order at f = 1, but of first 
order at small r. The FM-AM phase transition is always of first order. 

3.3 Perturbative determination of the counterterm coeffi- 
cients 

The counterterm coefficients k, Ai + A3 and A2 can be calculated order by order in 
perturbation theory by demanding the Slavnov-Taylor identity 

{A,ip) A,{q)) =^Sip + q) (3.4) 

to be satisfied in the continuum limit a — > 0. (The Slavnov-Taylor identity does 
not determine the other linear combination, Ai — A3, because it corresponds to a 
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Figure 2: T/ie phase diagrams at r = 1 (Figs, a and b) and r = 0.2 (Figs, 

c and d). The plots on the left were obtained by a mean-field calculation, and those 
on the right by Monte Carlo simulations on a 4^ lattice. The dash-dotted lines mark 
the symmetry line (cf. Eq. ( ^.18^ )). The perturbative result for the FM-FMD transi- 
tions is represented in the plots on the right by solid lines. The relation between the 



perturbative and Monte Carlo results will be discussed in Sect. 4-1 
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Note that the symmetry line in Figs, c and d coincides with the k, = axis. 
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Figure 4: The (k, f) phase diagrams for k, = 0.05 (Figs, a and b) and 0.2 (Figs, 
c and d). The mean-field phase diagrams are displayed again on the left, and those 
from the Monte Carlo simulations on the right. The dash-dotted line is the symmetry 
line. The perturbative result for the FM-FMD transitions is represented in the plots 
on the right by solid lines. The relation between the perturbative and Monte Carlo 



results will be discussed in Sect. 
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Figure 5: Same as Fig. but for k = 0.8. 



gauge-invariant operator.) 

First, we transcribe the left-hand side of eq. ( p.4|) to a finite lattice, 



where = 2sin(p^/2) is the lattice momentum, and 

\ 

. x,y 



W{p)=J2p,P^^Up)^ (3-5) 



Ajl(p) = ^2 13 rp felm Im [/^j^exp(ip {x - y))\ . (3.6) 



is the vector two-point function in momentum space, with L^T the lattice volume of a 
cylindrical lattice of spatial extent L and temporal extent T. The vector propagator 
A^^(p) is now computed to a given order in g'^. The resulting expression for A^^^{p), 
and hence also for W{p), will be a function of the counterterm coefficients. The 
counterterm coefficients are then determined such that in the continuum limit 

\imW{p)=^ (3.7) 

(keeping the physical volume fixed). Similarly, the coefficients A4 and A5 of the quartic 
counterterms can be calculated by requiring the Slavnov-Taylor identity 



= {5{P + q) S{k + l) + 6{p + k) 6{q + l) + 6{p + I) 6{q + k)) (3.8) 
to be satisfied in the continuum limit a — 0. 
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AM 



Figure 6: Schematic three-dimensional plot of the (k, k, f) phase diagram in the 
range from r = Otor = l,K = OtoK = l. The two dashed lines on the front faces of 
the phase diagram cube indicate the position of the symmetry surface (Eq. \2.18^ ) on 
those faces. The thick line represents the FM-FMD-PM tricritical line. 



After inserting 



A^,co = JI^xp (i p ■ x) exp {i p^/2) A^{p) 



(3.9) 



with J2p = (L'^T)'^ Yl^p, into the action (|2.3| ), we obtain from the term bihnear in 
vector potential Af^{p) for the tree- level vector propagator 

-1 



if + m^) V - f 1 - i j p^ Pu 



(3.10) 



where 



(3.11) 

We included the mass counterterm in the tree-level propagator, since it also functions 
as an infrared cutoff. 

We have calculated the critical values of the coefficients Ai, A2 and A3 to 
one-loop order in perturbation theory. To one-loop order the vector propagator is 



(m2 + (l + Ai/)p2 + A2/py 5, 
15 
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Figure 7: Projection of the FM-FMD-PM and FM-AM-PM tricritical lines onto a 
K = const plane obtained by a mean-field calculation (Fig. a) and by a Monte Carlo 
simulation (Fig. b). The FM-FMD-PM and FM-AM-PM tricritical lines in Fig. a 
appear to merge into an FM-FMD-AM tricritical line at large H. The two solid curves 
in Fig. b are to guide the eye. 
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v,(o) 



and 



k a 



(3.12) 
(3.13) 

(3.14) 



k a 



2 ^a/l°'' (^) ka 



k OL 



_9_ 
2 

-9' 



Y.^lf\k)-\7Y.^]^!P\k) k, k 

. k k 

HpIY. [kl A^\k) + kl ^T\k) - 2 ^J.Arf (A;) 

Pi^Pv ■ 



(3.15) 



(3.16) 



. k k 

^f^uip) T,^^{p) are, respectively, the contributions from the gauge- fixing and 
plaquette actions. The Z-factor, Eq. (|3.13| ), originates from the fact that we used the 
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composite operator Im U^x instead of A^^, in the definition of A^^{p). On a symmetric 
lattice {i.e. L = T) the self-energies are 

29' 



9' 



r/ii-i (/ii + 3/i2)+pJr [K 



^tip) = -r {J: 



/12 



12 



^11 

i f (111 - /12) 



(3.17) 
(3.18) 



(note that the last expression is transverse, as it should be) with lattice integrals 



k 



k 

The Slavnov-Taylor identity, Eq. { \i.7\ j, is satisfied to one-loop order for 

8 ' 



A. 
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and 



with 



X2 
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-FMD — — 






(e - 1)) ^2 


Ai + A3 = 


= -^i (1- 


— P 

k 1*1* " 


= 0.0951(1) , 




= 0.1549(1) , 


?2 t2 

\ ^ 1^1 1^2 


3 V4 



3 ; 



0.0516(1) 



(3.19) 
(3.20) 
(3.21) 

(3.22) 
(3.23) 

(3.24) 

(3.25) 
(3.26) 

(3.27) 



in the limits L, T — 00 and m? — > 0. Equation ( p.22| ) provides us with an expression 
for the FM-FMD phase boundary, to be compared (in Sect. |4.1|) with our Monte 
Carlo results. One can also verify that in the limit 9 ^ Q Eq. ( pj.22D turns into the 
corresponding result for the reduced model , taking the limit 9 ^ Q such that k in 
Eq. (|2.7|) is kept fixed, i.e. f — 00. (Note that both A2 and the sum Ai + A3 vanish 
in the limit r ^ 0. This happens because of a combination of two facts: for f = 
a ghost action can be added such that the full action has an exact BRST symmetry 
on the lattice however, to one loop, the ghosts do not appear in the vacuum 
polarization in the abelian case. The mass counterterm does not vanish for r = 
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Figure 8: The ratio W{p)/C, as a function of on a 6^24 lattice at the point 
{g,K,r) = (0.6,0.8,1). See text. 

since BRST invariance allows for equal non-zero masses for the U(l) gauge field and 
the Fadeev-Popov ghosts.) We can fix Ai and A3 if we demand the wavefunction 
renormalization constant to be equal to one at the one-loop level, then 

which then determines A3 using Eq. ( |3.24| ). 

As a check of our results, and to get a feeling for the effects of the various 
counterterms, we have calculated W{p) in Eq. ( |3.5| ) to one-loop order, and plotted in 
Fig. I the ratio W{p)/^ as a function of for < < vr, /i = 1, . . . , 4. After tuning 
the counterterm coefficients this ratio should approach one in the continuum limit 
a — >■ 0. The inversion of the 4x4 matrix in Eq. (|3.12|) was done numerically. The 
ratio W{p)/^ was computed on a 6^24 lattice at the point {g, K,r) = (0.6, 0.8, 1). We 
used antiperiodic boundary conditions in the time direction to avoid the zero mode 
of the propagator. The open triangles in Fig. ^ correspond to setting k = 0, Aj = 0, 
i = 1, . . . , 5. The plot shows that the ratio W{p)/C, is clearly below one and, as shown 
by the irregularities, is not a continuous function of p^. These irregularities are caused 
by the part of the self-energy with the structure of the A2 counterterm. Setting A2 
equal to the value obtained from (|3.23| ), one gets the values represented in Fig. || by 



the squares. All irregularities disappear and W{p)/C, is a continuous function of p^. 
The filled triangles in Fig. ^ were obtained after setting n = kfm-fmd (Eq. ( |3.22D ) 
and the crosses after setting also Ai + A3 to the value determined from Eq. ( |3.24| ). 
The graph shows the crosses to be very close to one for all momenta. 



3.4 Particle Spectrum at the FM-FMD Phase transition 

In this subsection we discuss the particle spectrum. We mentioned already in Sect. ^ 
that the particle spectrum in the FM phase of the pure gauge-Higgs model {R = 0) 
contains a massive vector boson and a massive Higgs particle whose masses scale when 
K is tuned towards the FM-PM phase transition. On the other hand, at the FM-FMD 
phase transition the action ( pl3|) is supposed to provide a new lattice discretization 
of a theory of free photons, with nothing else. A Higgs particle associated with 
the longitudinal gauge degrees of freedom should be absent. Therefore, the particle 
spectrum has to change when crossing the FM-FMD-PM tricritical line in the FM 
phase: an (unstable) Higgs bound state should exist near the PM-FM transition, and 
not near the FM-FMD transition. 

The question of the existence of a Higgs bound state is a non-perturbative issue. 
In this section, we present perturbative results for the various correlation functions 



that will be used to probe the spectrum numerically in Sect. ^]2. 

The vector propagator defined in Eq. ( p.6|) was already calculated to one-loop or- 
der in the last section, and is given by Eq. ( |3.12| ). It is evident that one indeed obtains 
a free, canonically normalized vector propagator if the four counterterm coefficients 
K, Ai, A2 and A3 are tuned towards the values given in Eqs. (|3.22| - |3.2^J3.28|) . 

An operator containing the quantum numbers of the Higgs particle is Re U^^- 
The corresponding Higgs two-point function on a cylindrical lattice is given by 

A^,{p) = j^(^Y.^eU,,ReU,yexp{tp{x-y))^ , (3.29) 
which has been used for the numerical determination of the Higgs mass in gauge-Higgs 



models [IS 



It is easy to verify that to one-loop order 

^fJ'\p) = \ 9' exp {^ {p, - p^)/2) £ Alf\k) A^f^p + k) . (3.30) 

k 

For small p, we can extract the non-analytic part by replacing the integrand with its 
continuum expression, and we obtain (for w? = and L,T ^ 00) 

EAJ^^Hp)! =-il9')j^,{S + e)\ogp\ (3.31) 

non— analytic ^ yATTj 

It is obviously not possible to conclude from this perturbative calculation alone that 
a Higgs particle does not exist. However, one may compare a nonperturbative evalu- 
ation of the same correlation function with the perturbative result. If they agree, this 
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provides evidence that a Higgs bound state does not occur in the theory, and this is 
what we will investigate in Sect. ^72 . 

Another equivalent way of looking at this is to consider the coordinate-space 
correlation function 

Cf.dx - 1/1) = (Re U^, Re U,y) - (Re U^,) (Re U,,) , (3.32) 

which, if no bound state is present in the spectrum, should factorize for |x — ?/| — > oo 

G;^Mx-y\) = C,, \GlA\x-y\)f , (3.33) 



as 



where 



GLi\x - y\) = (Im f/^x Im U^y) 



(3.34) 



is the vector correlation function, and Cfj,^ is a constant which can be determined in 
perturbation theory. 

To leading order in perturbation theory we find, 

4 4 
^^^u\\X ~ y\) — ~^\^iix ^uy/O — l^^^x.vy 

and 



(3.35) 



Gj'.dx - 1/1) = / (A^. A,y\ = g' A25, , (3.36) 

where (■ ■ •)o denotes the quantum average with the part of the lattice action (|2.3| ) 
that is quadratic in A^, and 



(3.37) 



with A^;,(o)(p) given in Eq. (|3lO| ). Substituting Eqs. and ( WMl into ( F33D 

leads to 

C^. = ^. (3.38) 



1 
2 



The expressions in Eqs. ( |3.35| ) and ( |3.36| ) are represented by Feynman diagrams la 
and 2a in Fig. respectively. Note that, to leading order in perturbation theory, 
factorization holds without any tuning of the six counterterm coefficients k, Ai, . . . , A5. 

We now wish to verify explicitly that factorization holds also to next-to-leading 
order in g^. To this order we obtain. 



G^i,{\x ~ y\) = ~^{Afix\yi^ " ))o — -^{Afix^uy + '^^xAuy)o 



9_ I ^y,(o) 2 [ ^ 



9^ 



fj.x,uy 



fix, fix 



uy,uy 



fix uy 
1 

2 



(^.^^J^Oo ,(3.39) 



where designates all terms of the lattice action ( |2.3| ) which are quartic in the 
gauge potential A^^^. Similarly, we find for G^j^{\x — y\) 



GU\^ -y\) = g'{A,.^A,y (1 - 5f )))o - ^{A,.^A% + A%A^y)o 



/\V,{0) 
fix,uy 



6 

A^.(O) _lA^.(0)' 
fix, fix ' uy,uy 



{A^^A^y 5;'))o ) . (3.40) 
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The various diagrams that contribute to and are displayed in Fig. |^. Diagrams 
lb, Ic, 2b and 2c correspond to the terms in Eqs. p.39| ) and ( ^.401 ) which are 
proportional to [A^^(°)^ + Ajf^(°y, and give only a contribution to the wave-function 
renormalization constant. The four-point vertices in diagrams Id-lf (we will refer 
to diagram If as the "figure-eight diagram") and diagram 2d arise from Si'^\ The 
integral expressions for those diagrams are given in Appendix B. 

In perturbation theory, one expects that — y|) ~ {G^j^{\x — y\)Y for large 

|x — Here, we show this to be true also at two loops. It can easily be verified that, 
after squaring G^^,{\x — y\), diagrams 2a, 2b and 2c combine into lb and Ic and that 
similarly 2a and 2d combine into Id and le. The Higgs correlation function can then 
be written as 

G'2.(|x-i/|)= fU/C«) \Gl{\x-y\)f-4{AlAlsi'yj + 0{g'), (3.41) 



= 1 (K,, + Ku) , (3.42) 



where 



[cf. Eq. (|]TD) and {Al^Al si^'^)l^ is the contribution from the fi gure-eight diagram 



which is given by the complicated expression in Eq. ( |B.4| ) in Appendix B. The various 
terms contributing to the figure-eight diagram in Eq. (|B.4]) can be divided into two 
classes according to the number of loop momentum factors in the loop integrals. First, 
each of the terms in Eq. ( p.4|) consists of the product of two one-loop integrals. The 
terms in Eq. ( |B.5D contain in one of the two integrands no explicit momentum factors 
{p etc.), whereas the terms in Eq. (|B.6|) contain in each of the integrands at least 
one such momentum factor. The momentum factors correspond in coordinate space 
to derivatives that act on that particular loop. A dimensional analysis then shows 
that the terms in Eq. ( [B.6| ) do not give a contribution at large distances (they give 
only contact terms in momentum space). The only contribution at large separations 
comes from the terms in Eq. (|B.5|) . The loop integrals without any explicit momentum 
factors behave at large separations |a; — 1/| as the leading order term {A'j^^Aly)^ which, 
in the limit m — 0, is logarithmic divergent in momentum space. In each term, the 
second integral contains two momenta and approaches a constant in limit k ^ 0, thus 
leading to a contribution to the constant C^^^. Therefore, for |x — ?/| — oo, 

- '-^{^Ksiyj = / C^:J \gU\^ -y\)Y + 0{g') . (3.43) 



As a check (for the case jj, = u), and in order to determine the constant C^J^, we have 
numerically computed the three ratios 

/a2 a2 c'(4)uf 



I a2 a2 c.(4)af,I 

r\\^-y\) = ^l T,o)2 > (3-45) 



fix fiy 1 /O 
g L-l^x,fiy 

12 a2 c(4)\if,n 



/ A2 A2 CW\J 

r'\\x-y\) = , (3.46) 

g L^fix,fiy 

21 





2.) G"^ 



(a) n (b) (c) 



(d) 



Figure 9: Feynman diagrams for the coordinate-space correlation functions G^^^dx — 
y\) andG]:Mx-y\). 

:{4)\lf / a2 a2 c{4)\lf,I „j / a2 a2 Q(4)\lf,II 



where (4.'4.^r)o and {A^^Al^Snr are given in Eqs. O 



]5D and ( p 6|) of Appendix B respectively. We have chosen x and y as two on-axis 
points. The ratios were computed on a symmetric lattice of size at the same point 
in the phase diagram {{g, k, K,f) = (0.4, 0.1, 0.2, 1)) where the numerical simulations 
were performed (c/. Sect. ^]2]^.) If factorization holds the two ratios r{\x — y\) and 
r^{\x — y\) should exhibit a plateau at large separations \x — y\. As an example we 
have plotted r{\x — y\) in Fig. [lO| a as a function of |a; — |/|/L for L = 14 (crosses) 
and 16 (triangles). The plot shows that we obtain indeed a plateau at large \x — y\. 
In Fig. p!0| b we have plotted the mid points, i.e. the values of r(L/2), r^(L/2) and 
r^^(L/2), as a function of 1/L^ = (a/LP^^*^ )^ where L^^^^- = aL is a fixed physical 
scale and a is the lattice spacing. We see that r^^(L/2) indeed vanishes in the limit 
a —>■ 0. In contrast, the ratios r(L/2) and r^(L/2) approach a non-zero value in this 
limit. The constant C^^ is then given by 



C(2 = i limr(L/2) = 0.12(1). (3.47) 



L 



For the constants C^j] and C^^ we find the values (for g = OA) 

C«= 0.3277(1) and C^^ = ^ + + /C^J = 0.5716(16) . (3.48) 

The above arguments do not lead to a constraint on the counterterm coefficients Ai, 
A2, A3 and K at this order. It is however clear that our above arguments are true only 
if A4 = A5 = 0, which is consistent to this order in g"^. At higher orders, factorization 
holds only when the counterterms A4 and A5 are tuned appropriately. (Note that then 
for any Ai, A2 and A3 the theory is free in the a ^ limit.) 
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Figure 10: Fig. a shows the ratio r{\x — y\) as a function of \x — y\/L for L = 14 
(crosses) and 16 (triangles), with periodic boundary conditions. In Fig. b the ratios 
r{L/2), r\L/2) and r"(L/2) are displayed as a function ofl/L^. The straight lines 
in Fig. b are to guide the eye. 



4 Numerical results 

In all our Monte Carlo simulations we set Ai = A2 = A3 = A4 = A5 = 0. The action 
depends then only on the four parameters g, k,k and f . We have seen in the previous 
section that, for the quantities we will consider, A4 = A5 = in perturbation theory to 
the order we have taken into account. This means that we can set them equal to zero 
also in our numerical computations, as long as the latter agree well with perturbation 
theory (within our precision) with the same choice of values. In addition, we will be 
mainly concerned with factorization (c/. previous section), which should work for any 
choice of Ai, A2 and A3. 



4.1 Phase diagram 

For the determination of the {k, k, f)-phase diagram we have to construct observables 
which allow us to locate the various phase transitions. The observables we used are 
three internal energies 

1 d 

En = T — In Z , (4.2) 
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g-f- 



7 



which in the vector picture are given by the expressions 



Ep 

Err.f. 



4L4 



\X,fJ.<U 



64L4 



(4.3) 

(4.4) 
(4.5) 
(4.6) 



These quantities are not order parameters, and hence do not vanish in any of the 
various phases, but they signal phase transitions by an abrupt change. In the case 
of a second order phase transition we expect to find, in the infinite- volume limit, 
an "S"-like curve with an infinite slope at the phase transition. At a first-order 
phase transition the internal energies exhibit a jump. On a finite lattice, however, it 
is difficult to distinguish between first- and second-order phase transitions, and it is 
usually necessary to perform a careful finite-size scaling analysis to settle the question 
of the order. In the FMD phase the hypercubic rotation invariance is broken by the 
non- vanishing vector condensate, (A^) 7^ 0. A true order parameter which allows to 
distinguish the FMD phase from the other phases can be defined on the lattice by 
the expression 



V 



A 



E 



(4.7) 



which reduces in the constant-field approximation to [| X]/^ sin^((?y4^)]^/^. On a small 
lattice, the system tunnels from one of the 16 discrete minima in Eq. ( p.3|) to the 
others. This is the reason why, in Eq. (|4.7|), we took the modulus of j^J^x^'"^ U^^- 
The summation over x is to project onto zero momentum. 

The Monte Carlo simulations were done with a standard 5-hit Metropolis algo- 
rithm, and were performed either in the vector or in the Higgs picture. We wrote 
two codes and checked that the results obtained in the two pictures are consistent. 
The vector-picture simulations require less CPU time since only the gauge fields have 
to be updated. However, the autocorrelation time for gauge non-invariant quantities 
turns out to be slightly larger for the vector-picture simulations. We furthermore 
find that the Higgs-picture simulations perform slightly better in regions of the phase 
diagram where metastabilities occur (the region near the FM-PM phase transition at 
f = and large k). Most of our simulations were carried out in the vector picture. 

We have explored the phase diagram again at g = 0.6, and, as in the mean-field 
analysis, we kept either the value of /5 or of f fixed and scanned the two-dimensional 
(k, f) or (k, k) plane. At each point of the scan we accumulated 5000 Metropo- 
lis sweeps, which were preceded by 3000 equilibration sweeps. The observables in 
Eqs. ( [4.4| - [4.7|) were measured after each sweep. We corrected for autocorrelation time 
effects by multiplying the statistical error bars with v^2r where r is an estimate for 
the integrated autocorrelation time, which for the local observables (4.4- 477|) is in the 
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Figure 11: Scans across the FM-FMD (Fig. a), FMD-PM (Fig. b) and FM-PM 

(Fig. a) phase transitions. The three curves in Fig. c which are labeled by the numbers 
(1), (2), (3) were obtained at {K,r) = (0.01,1), {K,r) = (0.05,0.15) and {K,r) = 
(0.05,-0.30). The estimated positions of the phase transitions are marked in all 
graphs by the vertical dotted lines. 
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range between 2 and 5. Most of the phase diagram scans were performed on a 4'^ 
lattice. A few runs at f = 1 were also done on an 8^ lattice. 



In Figs. |Tl|a-c we have displayed the various observables ( |4.4| - |4.7| ) for some exem- 
plary scans across the FM-FMD (Fig. |n]a), the FMD-PM (Fig. |n]b) and the FM-PM 
phase transition (Fig. pJ]c). Figure |ll]a shows that the order parameter V is very 
small in the FM phase and rises sharply at the FM-FMD transition. The fact that V 
is also non-zero in the FM phase is a finite-size effect. The internal energies and 
i?g.f. show a sharp kink at the phase transition, whereas Ep changes only very little. 
The transition seems to be of second order, in accordance with our mean field results 
and with perturbation theory. 

The position of a phase transition on a finite lattice can be defined in different 
ways [e.g. the position of the maximum of the specific heat, or the real part of the 
partition function zero with the smallest imaginary part), and these transition points 
will all differ slightly from each other by an amount which vanishes in the infinite- 
volume limit. In our case we have identified the FMD-FM phase transition from 
the point where the slope of V is maximal, but one should keep in mind that the 
position of the phase transition in the infinite-volume limit will slightly deviate from 
that value. 

The observables for a scan across the FMD-PM phase transition at f = 1 are 
shown in Fig. Ob. We believe that this phase transition is of second order which 



agrees with the findings from our mean-field calculation. The point where the slope 
of V is maximal always coincides nicely with the point where the slopes of E,^ and 
-Eg.f. are maximal. Finally, we have displayed in Fig. 11c the internal energy E^^ 



for three scans across the FM-PM transition ((1): f = 1, k = 0.01; (2): r = 0.15, 
K = 0.05; (3): f = —0.3, k, = 0.05) as a function of k. The plot indicates that the 
order of the FM-PM transition changes from second to first order when r is lowered 
from 1 to —0.3 at k > 0. This agrees with our mean-field calculation. In contrast 
to the mean-field calculation, our Monte Carlo simulations seem to indicate that the 
FM-FMD transition is of second order at small r. 

We have compiled our results for the various phase transitions into (k, /5)-or [k, r)- 
phase-diagram plots, which are displayed on the right in Figs. 0-^. We have again 
determined the phase diagram only above the symmetry lines (dash-dotted lines in 
Figs. 0-^), but we checked with a few scans that the phase diagram is indeed symmetric 
with respect to those lines. A schematic three-dimensional phase diagrams of the 
f > region is shown in Fig. Qualitatively, the Monte Carlo phase diagrams 
comply nicely with the mean-field results, which are displayed on the left in Figs. 



and which were discussed already in Sect, 



As mentioned already in Sect. |3.2| , a difference between the mean-field and Monte 
Carlo phase diagrams is observed at f = (Figs. |^ and d). The mean-field approxi- 
mation predicts an FMD phase at /t^O.25, |fi:|^0.125. This is not confirmed by the 
Monte Carlo simulations which give strong evidence that this region is filled by a PM 
phase. The Monte Carlo simulations indicate that the PM phase extends to very large 
values of R but, on the basis of the simulations, it is of course impossible to decide if it 
ends at a finite value of k, or if it extends to k = cxd. Figure ^ shows that the observed 
difference at f = is connected to the fact that in case of the mean-field calculation 
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the FM-FMD-PM tricritical line (in Fig. |^ we displayed a projection of that line to a 
constant-K plane) penetrates through the r = plane and at k ^ 0.3 merges with the 
FM-AM-PM line into a FM-FMD-AM line, whose projection approaches asymptot- 
ically the f = axis at large k. The FMD phase extends therefore slightly into the 
f < half space. In contrast, in the case of the Monte Carlo simulation, we find that 
the FM-FMD-PM line stays in the f > half space, and approaches the f = plane 
from above. The FMD phase resides only in the r > half space. The solid curves 
in Fig. 03 are to guide the eye, and were obtained by fitting the empirical ansatz 
r = ±a/{li + b) to the data with a and b constants. It turns out that the constants a 
and b are very similar for the FM-FMD-PM and FM-AM-PM lines. 

The continuum limit relevant for the gauge-fixing approach should be performed 
by approaching the FM-FMD phase transition from the FM side (away from the 
tricritical line at which the FM-FMD transition surface ends, for instance at r 1). 
Our Monte Carlo simulations indicate that this transition is of second order. In the 
next section, we will show that the particle spectrum in this continuum limit contains 
only a massless vector particle, the photon. The FM-PM phase transition at small 
r near the FM-FMD-PM line and the whole FM-AM phase transition seem to be of 
first order, and no continuum theory can be defined at those transitions. 

Finally, we compare our one-loop result for kfm-fmd (Eq- (|3.22| )) with our Monte 



Carlo results. The one-loop results for kfm-fmd are represented in the graphs on the 
right in Figs. by the solid lines. They agree reasonably well with our Monte Carlo 
data. In the region of the FM-FMD phase transition we find that the departure of 
the perturbative results from the Monte Carlo data is in all cases smaller than two 
standard deviations. In some cases (see for instance the phase diagram at f = 1 in 
Fig. Pd), however, the one-loop curve is systematically above the numerical results. 
To see if this deviation is due to finite-size effects we repeated some of the scans across 
the FM-FMD phase transition on an 8^ lattice. The obtained transition points are 
marked by the open circles in Fig. ^d. They do not significantly deviate from our 
data on the 4^ lattice. We also evaluated the lattice integral Xi in Eq. ( |3.22| ) on a 
4^ and 8^ lattice which again did not result in a significant change of the situation 
encountered in Fig. ^. We therefore believe that the observed deviations are due to 
higher loop corrections. They become smaller for smaller r. This is not unreasonable: 
as long as perturbation theory applies, some couplings are proportional to f, so that 
higher loop terms maybe less important for smaller f . (Of course, when f is too small, 
perturbation theory breaks down all together, as is clear from the case f = 0, where 
we do not even find an FM-FMD transition, numerically.) 

4.2 Vector and Higgs two-point functions 

To see whether the spectrum at the FM-FMD phase transition indeed contains only 
a massless photon we have computed the vector and Higgs two-point functions nu- 
merically. 
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Figure 12: The inverse vector two-point function (points with error bars). The long- 
dashed line is the tree-level result, and the solid line the one-loop result. The latter 
is almost on top of a linear fit (short dashes) to the data. Fig. a shows results for 
K = 0.01, while Fig. b shows a blow-up of the small momentum region for k =0.8, 0.4 
and 0.01. 



4.2.1 Vector two-point function 



In our Monte Carlo simulations we computed the correlation function A^^(p) in 



Eq. ( |3.6| ). We have set = z/ and p^ = 0. The two-point functions A^^{p), = 1, 2, 3 
were determined for a table of lattice momenta p which lead to different values of 
p^ = J^iy^fiPl- Finally, we took the average of the three two-point functions A^^(p), 
H = 1,2, 3. All Monte Carlo simulations were performed at Aj = 0, z = 1, . . . , 5. 

For the above choices, we first derive a simplified expression for the vector two- 
point function to one-loop order, which can then easily be compared with the Monte 
Carlo results. After setting /i = i/ = 1 and pi = 0, the one-loop vector two-point 
function in Eq. ( |3.12| ) can be written as. 



Ari^'\p) 



A^'(o)(p)-i + S(p) 



11, Pi=0 



^ii(/) fA^'(°)(p) - A^'(°)(p) A^'W(p) + 



11, pi=0 



with 



(4.8) 
(4.9) 

where Zii{g^) is the wave- function renormalization constant in Eq. ( |3.13| ) and Ti^^'^'^{p) 
is the self-energy in Eq. (|3.14| ) (m^ = 2g'^K is already included in A^'(°), cf Eq. (CT) )- 
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After writing Eq. as J^p^p) = a^{p)5^u+h^,{p)p^p, and Eq. (|3ig) as Aj^(o)(p) = 

•A.{p)5fj_y + B{p)pfj,pu, and using the fact that pi = 0, we obtain 

^ii^'\p) = ^ii(/)[^(p) + ^(p)«i(p)^(p) + ---U=o 

= . (4.10) 

Note that the expression in Eq. ( ^.101 ) does not depend on A2 and A3 as we set pi = 0. 
This also imphes that we will not see any effects of the Lorentz-symmetry breaking 
part of the one-loop self-energy in our data. In the following we set also Ai equal to 
zero since in our Monte Carlo simulations we have also not included this counterterm. 
On a lattice which is asymmetric in space and time ai{p) is given by the expression 

2 ~ 2 
«i(p)U=o = ^ [3Ai + lu] - ^ [2/12 + /l4 + hi] 



-| (7' [2 (J12 - /12) (pI+pD+pI {Ju + J41 - 2 Jm))] , (4.11) 

where the integrals /^i/ and J^i, are given in Eqs. (|3.2CI|) and ( p.21|) . 

In Fig. |T2| we plotted the results of our numerical computation of the inverse 
vector propagator in momentum space, as a function of The long-dashed and 
solid lines represent tree-level and one-loop perturbation theory evaluations of the 
same quantity, at values of the parameters equal to those used in the numerical 
computation. In Fig. [T^, these are (K,K,f) = (0.01,0.2,1), while in Fig. \r2j o we 
show an enlargement of the small-momentum region for {K,r) = (0.2,1) and three 
different values of k, 0.8, 0.4 and 0.01, the latter being very close to the FM-FMD 
transition point. A linear fit (short dashes) to the numerical results can, within 
the resolution of the plots, hardly be distinguished from the one-loop curve. The 
gauge coupling g is 0.4 and the lattice size is 6^24. The perturbative expressions for 
(the inverses of) A^'^^\p) and A^'^^\p) were evaluated on a lattice of the same size. 
We have measured the vector two-point function A^^(p) on 2 x 10^ configurations 
which were generated by our 5-hit Metropolis program. The error bars were again 
computed by multiplying the standard deviation with ^/2r where r is an estimate of 
the integrated autocorrelation time obtained from the autocorrelation function. 

From these results we draw two conclusions: the fact that a linear fit works very 
well confirms that the theory is a theory of free photons near the FM-FMD transition, 
and the good agreement with perturbation theory implies that this can be understood 
in perturbation theory, as explained in Sect. |^. 

For comparison, we have repeated this analysis at a series of points near the FM- 
PM phase transition, at k > and f = 1. We found that the perturbative results 
do not converge (one loop is not close to tree level), and also do not describe the 
numerical data in this case (with deviations well over 100%). We also determined the 
vector boson mass my by fitting an ansatz Zyjijp' + m^) to the numerical data of 
the vector two-point function at small We find that the obtained vector boson 
mass my shows qualitatively the same k dependence as reported in Ref. [|1^] for the 
U(l) gauge-Higgs model. The vector boson mass decreases when k is lowered from 
the FM side towards the FM-PM transition. However, it does not vanish at the phase 
transition which is probably, like in the U(l) gauge-Higgs model, a finite-size effect. 
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Figure 13: The inverse Higgs two-point function, A^(p) , plotted as a function of 
IpP for a series of points near the FM-FMD phase transition at g = 0.4, k = 0.2 and 
r = 1. The error bars mark the Monte Carlo results. The lattice size is again 6^24. 
The perturbative results (Eq. (\3. 3(\ ) ) are represented by the solid lines. 



4.2.2 Higgs two-point function 

As mentioned before, the spectrum at the FM-FMD phase transition should contain 
only a massless photon, and no Higgs particle. In this section we present the Monte 
Carlo results for the Higgs two-point function, and compare them with the analytic 



formulas derived in Sect. 3.4 



We have computed the momentum space Higgs two-point function in Eq. ( 3.29| ) 



As in the case of the vector two-point function, we have set z/ = /i and p^ = 0. The 
Higgs correlation function was determined for the same lattice-momenta as the vector 
two-point function. 

To see if the two-point function leads to a pole we have plotted in Fig. A^(p) 
as a function of |pp = ^^^^pl for several k values near the FM-FMD phase transition. 
The simulations were performed a.i g = 0.4, and at k = 0.2 and f = 1. The numerical 
data are represented in Fig. |T^ by the error bars. The five data sets from the top to 
bottom correspond to k = 0.8, 0.4, 0.3, 0.1 and 0.01. At each k we have measured 
A^(p) on 2 X 10^ equilibrium configurations which were again generated with the 
5-hit Monte Carlo algorithm. As in the case of the vector two-point function, we 
corrected for the autocorrelation time effects by multiplying the standard deviation 
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Figure 14: The ratio R{\x — ?/|) as a function of \x — y\ at a point near the FM- 
FMD phase transition ({g,K,,K,,r) = (0.4, 0.1, 0.1, 1)^). The lattice size is 6^24. The 
horizontal dashed and solid lines were obtained by evaluating the ratio in Eq. ( H-lSj ) 
in the infinite volume limit to leading and next-to-leading order in perturbation theory 
(Eqs. (\3. 3^ and ( ^.48^ )-) The triangles were obtained by evaluating G^^{\x — y\) and 
Gj^^{\x — y\) in Eq. to next-to-leading order on the same lattice as used in the 

simulations. 



with a factor v2r. Our Monte Carlo simulations indicate that the FM-FMD phase 
transition is situated at k ^ 0. 

If the pole scenario would be correct, the A^(p) data at sufficiently small 
momenta should fall on a straight line, for k \ kfm-fmd- What we find, however, 
is that, when lowering k towards the FM-FMD phase transition, a cusp emerges at 
small momenta. Evidently, the data at small do not fall on a straight line. The 
cusp is due to the logarithm in Eq. ( p.31 ). The solid lines in Fig. 13 were obtained by 
evaluating Eq. ( p.30| ) on a lattice of the same size and for the same parameter values 
as used in the simulations. We find that the perturbative formula (|3.30|) describes 
the data very well. (A similar behavior was discovered before in the reduced limit 
of the model at f = 1 for the left-handed neutral and right-handed charged fermion 
propagators, which also exhibit such a logarithmic singularity at small momenta, and 
do not exist as bound states, see Refs. 



10,0 



Again for comparison, we looked also at the Higgs two-point function near the 
FM-PM phase transition at k = and g = 0.6. We find that kfm-pm ~ 0.18. We 
find that in this case, in accordance with expectations, the spectrum at the FM-PM 
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Figure 15: The ratio R{\x — y\) as a function of \x — y\ at a point near the FM-PM 
phase transition ({g,K,,K,,r) = (0.4, 0.08, 0.15, The lattice size is 6^24. The hori- 
zontal solid line and the triangles were obtained by evaluating the ratio in Eq. ( H-lSj ) 
to leading and next-to-leading order in perturbation theory. 



phase transition contains a massive Higgs particle, giving rise to a pole in A^(p). 
To extract the Higgs mass we fitted the A^{p)~^ data at small momenta to the 
ansatz {rajj- + \p\'^)/Zh- The Higgs mass decreases when k is lowered toward the 
FM-PM phase transition but, because of finite size effects, does not vanish at the 



phase transition (again, as in Ref. |T^). We also find that, as in the case of the 



vector two-point function, perturbation theory does not describe the PM side of the 
FM-FMD-PM tricritical line. 



4.2.3 Factorization 

The results of the previous section already give strong evidence that the Higgs two- 
point function in the continuum limit at the FM-FMD phase transition does not have 
a pole, and that consequently a Higgs particle does not exist in the spectrum. In 
Sect. p.4| we have shown that, to next-to-leading order in perturbation theory, the 
Higgs two-point function in coordinate space factorizes into the product of two vector 
two-point functions (Eq. ( |3.33D ). For A4 = A5 = 0, factorization holds to this order 
for arbitrary values of the other counterterm coefficients. 

In this section we will investigate whether factorization also holds non-perturba- 
tively, which would imply that a Higgs bound state does not exist. To this end, we 
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have computed the Higgs and vector two-point correlation functions G^^{\x — y\) and 
G^^{\x — y\) in Eqs. ( |3.32| ) and ( 13.34] ) in our Monte Carlo simulations as function of 



|a; — where x and y were chosen to be two on- axis points. The simulations were 
carried out at the point {g,K,R,r) = (0.4,0.1,0.2,1), which is located in the FM 
phase near the FM-FMD phase transition. As in our other Monte Carlo simulations, 
we have set Ai = . . . = A5 = 0. (Strictly speaking, we expect factorization only to 
hold for the properly tuned values of A4 and A5, but to the extent that our simulation 
results agree with our perturbative results, we do not expect to see the difference, cf. 
end of Sect. |^.) The lattice size is again 6^24. We find that enormous statistics is 
required to obtain a stable signal for G^^{\x—y\). This is partly due to the subtraction 
in Eq. (|3.32|). We used translation invariance on the lattice to improve the signal. 



We have accumulated in total about 27 x 10 Metropolis sweeps. The Higgs- and 
vector-correlation functions were measured after each sweep. 



In Fig. 14 we plotted the ratio 



G^,{\x-y\) 
G]:,{\x-y\) 



R{\x-y\) = -^ ^ (4.12) 



as a function of \x — y\. The error bars of the ratio were calculated by a blocking 
procedure. Figure |T^ shows that the Monte Carlo data for R{\x — y\) (crosses) are, 
within error bars, independent of \x — y\, indicating that factorization holds also non- 
perturbatively. In fact, factorization sets in almost immediately (both perturbatively 
and numerically). 



From our calculation in Sect. |3.4| it follows that this ratio is equal to 1/2 to leading 
order in perturbation theory. Figure |1^ shows that the Monte Carlo results are indeed 
very close to this value, marked by the dashed line. We find however that the data 
at small distances where the error bars are smallest are systematically above this 
value by a small amount. To understand this small discrepancy, we have computed 
the ratio also to next-to-leading order in perturbation theory. The next-to-leading 
order result for C^^, which was computed in Sect. |3.4| in the infinite volume limit, is 
represented in Fig. |n| by the horizontal solid line, which is indeed much closer to the 
numerical data. In addition, we have also evaluated all Feynman diagrams in Fig. ^ 
numerically on the same lattice and for the same parameter values as used in the 



Monte Carlo simulations. Figure |n| shows that the next-to-leading order result for 
R{\x — y\) (triangles) agree within two standard deviations for all separations |x — y| 
with the Monte Carlo data. We consider the \x — ?/|-independence of our Monte 
Carlo results for R{\x — y\) and the good agreement with the perturbative result as 
strong indication that factorization holds, and that consequently the spectrum near 
the FM-FMD phase transition does not contain a Higgs particle. We repeated the 
same calculation at a point in the FM phase near the FM-PM phase transition where 
a Higgs particle is known to exist. The Monte Carlo simulations were performed at 
{g, K, K, f) = (0.4, 0.08, 0.015, 1) with less statistics, about 5 x 10^ Metropolis sweeps. 



Figure [T5| shows clearly that, in contrast to Fig. |1^, the ratio R{\x — y\) depends 
strongly on |x — ?/| (notice also the difference in ordinate scale between Figs, [l^ and 
|T3|!), and seems to oscillate when the separation |x — y| is increased. The solid lines 
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and the triangles represent again the leading and next-to-leading order perturbative 
results for the ratio, both far off from the Monte Carlo data. It is obvious that 
factorization does not hold on the PM side of the FM-FMD-PM tricritical line. 

5 Conclusion and Outlook 

In this paper we investigated the gauge sector of the gauge-fixing approach for the case 
of a U(l) gauge group. This approach provides a completely new non-perturbative 
formulation of a U(l) gauge theory on the lattice, which is more complicated than 
Wilson's manifestly gauge-invariant compact formulation, but closer in spirit to the 
continuum formulation, for which gauge-fixing is indispensable. Gauge fixing allows 
us to control the longitudinal degrees of freedom which otherwise, as we mentioned in 
Sect. D, form a central obstruction to the construction of lattice chiral gauge theories. 

The action is rather complicated in comparison with the Wilson plaquette action. 
Apart from the Wilson plaquette term and the gauge-fixing term, it also includes 
six counterterms whose coefficients have to be adjusted such that Slavnov-Taylor 
identities are restored in the continuum limit. We have demonstrated in this paper 
that this new lattice formulation reproduces the desired properties of the continuum 
theory at a new type of continuous phase transition, the FM-FMD transition. The 
spectrum contains at this phase transition only a massless photon, while a Higgs-like 
excitation, associated with propagating longitudinal gauge degrees of freedom, does 
not exist. 

This new phase transition occurs in a region of the phase diagram accessible 
to weak-coupling lattice perturbation theory, as one would expect from the close 
relation to the continuum theory. In perturbation theory, by construction, the theory 
(without fermions) at the phase transition is a theory of free massless photons. The 
very good agreement between one-loop perturbation theory and numerical results 
makes us confident that this is also true non-perturbatively. 

As mentioned above, in order to make this work, six counterterms need to be ad- 
justed. However, only one of those (a gauge-field mass term) has dimension smaller 
than four. One expects that the others can be reliably calculated in perturbation the- 
ory. Our results indicate that this is indeed true: within our numerical precision, we 
find that we can do with just the tree-level values of all dimension four counterterms. 
Adding fermions to the theory will require a few additional counterterms, but (using 
a fermion formulation with shift symmetry [^) all of those are dimension four, and 
we expect that they can be calculated in perturbation theory as well. This means 
that the fact that counterterms are required does not make this formulation of lattice 
chiral gauge theories particularly expensive. 

It is clear that the precise form of the lattice action should be chosen such that 
lattice perturbation theory applies. It is therefore important to construct the lattice 
gauge-fixing action such that the dense set of lattice Gribov copies of the perturbative 
vacuum (U^^ = 1), which occurs for a naive discretization of the gauge-fixing action, 
is removed by adding higher- dimensional operators (with dimension six or higher). 
For comparison, we studied also the limit where those higher- dimensional terms are 
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omitted (by setting f = 0), such that these lattice Gribov copies are present. Our 
numerical results for this naive choice show that, at small k, there is an FM-PM type 
phase transition, in a universality class different from the FM-FMD transition (there 
is even some evidence that the FM-PM phase transition is of first order, implying 
that a continuum limit cannot be performed at all). The situation at large k is very 
unclear since the FM-PM phase transition is "wedged" (in the {k, k, r) phase diagram) 
between two tricritical lines, resulting in a very complicated phase structure, where 
four phases get very close to each other. We furthermore find that the numerical 
simulations in that region are hampered by strong metastabilities. These findings 
strongly suggest that for r = [i.e. naive gauge fixing) no phase transition in the 
desired universality class of the continuous FM-FMD transition (found at f ^ 1 and 
large enough k) occurs. This implies that the naive, f = action does not lead 
to a theory of free photons, and is unsuitable for the construction of lattice chiral 
gauge theories. It is likely that this is related to the dense set of lattice Gribov copies 
which occurs at r = 0, since they represent unsuppressed rough fluctuations of the 
longitudinal gauge field. In addition, our mean-field results indicate that for small 
nonzero r the FM-FMD transition may become first order, which would imply that 
small values of f should be avoided altogether. 

Our previous results on the fermion spectrum in the reduced model [l^ , combined 
with the results of this paper, provide what we consider to be convincing evidence 
that the gauge-fixing approach does indeed lead to a viable non-perturbative lattice 
formulation of chiral gauge theories for abelian gauge groups. We would like to 
emphasize that a key element of this approach lies in the fact that lattice perturbation 
theory provides a valid approximation to our lattice theory (including fermions 

0)- 

We believe that, in addition, our gauge-fixed lattice formulation has matured to 
the point where it may be used as an alternative to other gauge-fixing methods for 
abelian theories. In the traditional approach, one first performs a Monte-Carlo update 
using only the gauge action (Eq. (P^)). Then, a sequence of gauge transformations is 
performed, aiming to find "the best" configuration on the same orbit. For example, 
in the Landau-gauge method one attempts to maximize Re X]a;,At f^.M- ^ well-known 
problem is that local algorithms cannot (and do not) always find the global maximum. 
Some specific obstructions of global nature have been described in the literature (see 
e.g. Ref. and references therein). 

In contrast, in our approach one always performs the Monte-Carlo update on the 
entire gauge-field space, with one single Boltzmann weight. While the action ( p.3| ) is 
more complicated, this approach may nevertheless have several advantages: 1) one 
has all systematic errors completely under control; 2) we believe that global features 
{e.g. Double Dirac Sheets ||22|) do not cause any special difficulty in our approach. 



and that this is related to the fact that at no stage is our Monte-Carlo update con- 
strained to stay only on a single orbit; 3) finally, as mentioned earlier, we find that 
our photon propagator is in excellent agreement with theoretical expectation. (In 
order to maintain the good agreement of the photon propagator with the continuum 
theory for all four- momenta we expect that, of the five marginal counterterms, at 
least the one-loop A2 counterterm should be included in the Monte-Carlo update, see 
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Sect. ^73| , in particular the discussion of Fig. |.) 

Coming back to the program of constructing chiral lattice gauge theories, as a next 
step the gauge-fixing approach should be generalized to non-abelian gauge groups. 
This, even without fermions, is a difficult task. The main obstacle is the existence of 



"continuum" Gribov copies [^, which occur in non-abelian theories, in addition to 
the lattice Gribov copies discussed in this paper. It is well-known that the determinant 
of the Fadeev-Popov operator can be negative for some of these continuum Gribov 
copies, giving rise to a non-positive integration measure. It is therefore very likely 
that the weighting of gauge configurations in the path integral will deteriorate due to 
the presence of these continuum Gribov copies. In the worst case, contributions from 
different Gribov copies may even cancel each other (this is in fact not unlikely, in 
view of Neuberger's theorem pOf). This obstacle may be circumvented by adopting a 
gauge-fixing procedure as proposed in Refs. El, EBI. An advantage of this procedure 



is that the (gauge-field) integration measure is guaranteed to be positive. A possible 
disadvantage of this method is the fact that the counterpart of the Fadeev-Popov 
action is a highly non-local functional of the gauge fields. This makes a perturbative 
analysis, and, in particular, the construction of the counterterm action non-trivial. 
Work on this is in progress. 

Another project for future investigation concerns fermion-number violation. Most 
lattice chiral fermion actions (including that of Ref. [1^) can be written in the form 
J2x,y'4'xDx,y{U)'il'y with D^^yiU) the lattice Dirac operator. Obviously, this action 
(and also the fermion measure) are invariant under an exact global U(l) symmetry 
which, at first glance, seems to be in contradiction with fermion-number violation 



However, Ref. demonstrated, in a two-dimensional toy model, that fermion- 
number violation can still occur despite this exact symmetry. The central observation 
is that fermionic states are excitations relative to the vacuum. The global U(l) 
symmetry prohibits a given state to change fermion number, but nothing prevents 
the ground state to change when an external field is applied. We expect that a similar 
phenomenon may explain how fermion-number violating processes take place in our 
four-dimensional dynamical theory. 
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Appendix A 



We look for a translation-invariant solution, choosing the mean-field ansatz 

(l>x = f , (A.l) 

t/^a; = u exp {i A^) . (A. 2) 

The corresponding magnetic fields were replaced by 

hx = h^ , (A. 3) 

and 

H^^ = K exp {i A^) . (A.4) 

The 4 + (i mean fields m, /i<^, hu and A^, /i = 1, . . . , c? are space-time independent. 
Using these expressions we obtained for the free energy 

J^(v?, h^, u, hu, A] K,K,r) = L'^ <2(ph^ + 2duhu- log h{2h^) - d log /o(2/i„) 



^2.^«^^^^.^^~ ^ 



where L is the extent of the rf-dimensional lattice in each direction, and 

A; f) = -2 (4 d K f + k) u F{A) + k (1 +f) {2F{Af - d) 

+ 1. ^fM^ F(2A) (2c/+ 1) , (A.6) 
lo 

f'^\u, A] K,K,¥) =2ku'^ F{Af {¥ -I) - [q u^F{2Af 

-4 (2 1) + (rf- 1) M^) F(2A) 

+4c/(l + 2c/) + c/(2rf-5)M^] , (A.7) 

A; /€, K,r) = -^Kf (2 F{2A)^ - 2 (d - 1) F{2A) - d) , (A.8) 
16 ^ ^ 

f'\u,A; K,K,r) = ^ ™' (i^(2A) - d)' . (A.9) 

The two quantities F{A) and /o(/i) are given by 

F{A) = J2cosA^ (A. 10) 



and 



1 /■'^ 

/o(/i) = — da exp(±/i cos a) . (A. 11) 

TT Jo 



It can be checked that the above expression for JF reduces in the limit u —>■ 1 and 
f — > 1 to the free energy which we obtained before in the reduced model [||. 

The values of the 4 + d mean fields which are realized at a given point in the {k, 
K, r) parameter space are determined from the absolute minimum of the free energy. 
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The local extrema of the free energy are obtained by solving the A + d saddle-point 
equations 



dhu 
dT 



LP - 



U 











dA 
dJ 

dp! 







L'^ {C{(p,u; K, K,r) + B{(p^u; K, K,r) cosA^} sinA^ 
du 







where 



h{h) 



dlojh) 
dh 



(A.12) 
(A.13) 
(A. 14) 
(A.15) 

(A. 16) 



and 



C{ip, u; K, K, r) 
B{ip, u] K, K, r) 



(A.17) 



-2 K M v?^ + 4 K (1 + f) F{A) 
+4 K (f - 1) p?^ F{A) -Sdnrup?'^ , 
+ i K F (1 + 2d) pi^ -k¥ [\ (3 F{2A) 
-(2 {d+l) + {d- 1) u^)) + \ (2 F{2A) -id- 1)) p>^ 
+\ (F(2A) - rf) M^s} . (A.18) 



Equations (|A.14| - |AT5D can have multiple solutions and it is therefore important to 
insert the various solutions back into the free energy to find out which of them cor- 
responds to the absolute minimum. Since it is difficult to find a closed expression for 
the solutions of the saddle-point equations ( [A.12[ - [ArT5D , we minimized the free energy 
numerically. We set d = A. 

This was done in the following steps: The two saddle point equations ( |A.12D and 
( |A.13| ) were used to express the magnetic fields and hu in terms of p) and u. The 
free energy depends then only on p>, u and A^. The saddle-point equation ( [A. 14] ) has 
the solutions 



A^ = {z, z, z, z), {±A, z,z,z),..., {±A, ±A, ±A, ±A) , (A.19) 

where z = 0,n and cos A = —C/B (provided \C / B\ < 1). We have inserted each of 
these solutions back into the free energy JF in Eq. ( |A.5| ), which is now only a function 
of ip and u. Note that in the PM phase, where p> = 0, A^^ remains undetermined 
because all the dependence of the free energy on A^ {cf. Eq. ( [A.5D ) disappears. The 
minimization with respect to p> and u was done numerically by discretizing the (p), u)- 
space by a fine grid and calculating the free energy on each site of that grid. Finally 
we have picked the absolute minimum among the various A^ solutions in Eq. ( |A.19|) . 
The only solutions for which lead to an absolute minimum of the free energy turn 
out to be A^ = (0, 0, 0, 0), (vr, vr, vr, vr) and {±A, ±A, ±A, ±A) (with A^0,7f), which 
correspond respectively to FM, AM and FMD phases. 
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After this procedure we obtain at a given point in the phase diagram a numerical 
value for each of the eight mean fields u, h^, and A^, = 1, . . . , 4. The phase 
transitions, finally, were located by monitoring u, ip and A as a function of the four 
coupling constants g, k, k, and r. For locating the FM-PM and FMD-PM boundaries, 
we used the fact that (f = on the PM side and (f > on the FM and FMD side 
of the transition. Similarly, the FM-FMD (AM-FMD) phase transition was located 
using the fact that A = {A = n) in the FM (AM) phase and A 7^ 0, vr in the FMD 
phase. 



Appendix B 

In this appendix we give the explicit expressions for (A^a;^/xy5'f^'')o and {A'jj^^A'jj^ySi'^^)o. 
We find 



{A,.A^yS^\ = g'Y. e^'^^"^^ E Kf^ i^) (A^)A^(°) (k) (B.l) 



per 



and 



K^ax^uy'^l 10 — K^nx^uy^l /O ^ K^ax^vy"^! /O ^ \^ax^uy'-'l /O 



(B.2) 



where {A%Al^st^)l\ {A%Al^s['^)l^ and {A%A%s['^)l' are the contributions that 
correspond to the Feynman diagrams Id, le and If in Fig. ^ 



ia2 a.2 c;.(4)\ld_/.2 .2 0(4)^6 _ 2 



(4)\le 



J k {x-y) i(fc,,-fc,y)/2 



E ^If^ (p + fc) X E ^If ip) Sp.(p)A^(°) ip) 



pa 



/A2 a2 oW\1{_Ia2 a2 c(4) \ If ,I , / /i 2 a2 oW\ 
K^ax^uy'^l lO — K^ax^uy^^l 10 + K^ax^uy'^l I 



Ia2 a2 c;-(4)Uf,I 
K-^iix^vy'^l 10 



J k [x-y) i {kf_,-kv)/2 



X 



(B.3) 
(B.4) 



P V 



xEEA:r(g)Arr(g+fc)cosf COS 



^E EA,^) (p)A;f {p + k)p,p. 



p\ y p 



Ia2 a2 cW\^ 
K^lix^vy^^l h 



J k (x-y) i{k^-k„)/2 



Y { ^ E E ip) ^If ip + k)p,[p + k 



p p 
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(B.5) 



A q 

+4 E E Kf^ ip) Kf (p + k) h (p^k) 

pX P 

+^ E E Kf' (p) Kf ip + k) Pp cos ^ 



pA P 



j:KfHl)Kfi<l + k)g,cos^-^^ 



^ E E Kf iP) iP + k) PX {P'^k 

^ pX P 



^Y.^lf i<i)^lf {q + k)q,(q + k 

+^ E E ip) Kf^ (p + k) PX {pTk) 

^ pX P 

^T.^lfi<l)^:fi<l + k)q,[qTk) 
q 

+ E E Kf^ ip) ip + k) PX (p^k) 

pX P 

xE^T (g)A^,(°) iq + k)q,{qTk)^ 
q 

-2 E E Kf' (p) (p + k) PX {pTk) 

px p 

x£ Alf (q) {q + k) q, {qTk) ]) + ( ^ u 



X 



(B.6) 



where the self-energy S^^'''^(p) in Eqs. (|B.1|) and ( [B.3|) is given in Eq. ( p.l4|) . The terms 



in Eqs. ( p.5|) and (p.6| ) which are proportional to 1/^ are the contribution from the 
gauge-fixing action and all other terms arise from the part of the plaquette action. 
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